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We present a theoretical framework for the calculation of charge transport through nanowire-based 
Schottky-barrier field-effect transistors tliat is conceptually simple but still captures the relevant 
physical mechanisms of the transport process. Our approach combines two approaches on different 
length scales: (1) the finite elements method is used to model realistic device geometries and to 
calculate the electrostatic potential across the Schottky-barrier by solving the Poisson equation, 
and (2) the Landauer approach combined with the method of non-equilibrium Green's functions 
is employed to calculate the charge transport through the device. Our model correctly reproduces 
typical I-V characteristics of field-effect transistors and the dependence of the saturated drain current 
on the gate field and the device geometry are in good agreement with experiments. Our approach 
is suitable for one-dimensional Schottky-barrier field-effect transistors of arbitrary device geometry 
and it is intended to be a simulation platform for the development of nanowire-based sensors. 



I. INTRODUCTION 

Over the last few decades one-dimensional (ID) semi- 
conducting silicon-nanowires (SiNWs) have been widely 
studied as potential building blocks for future electronic 
devices due to their excellent electrical performance, 
small size, and controllable bottom-up fabricationii"— 
Many types of SiNW devices such as ultrasensitive 
sensors^ photodetectorsj^ and bipolar field effect tran- 
sistors (FETs)^i^ have been reported. The key issue for 
sensor applications is the down-scaling of FETs to one 
dimensional structures, such as nanowires. The attrac- 
tive feature of nanowire-based FETs is that the bind- 
ing of charged species can be directly monitored by the 
change in current through the NWs because of their high 
surfacc-to-volume ratios and small cross-sectional con- 
duction pathways. 

Recently Weber et al. have reported dopant- 
free Schottky-barrier (SB) FETs consisting of intrin- 
sic SiNWs working as a channel and NiSi2 nanowires 
working as source and drain contacts with gate lengths 
down to subphotolithographic values >i2rJ^ Their devices 
are based on single-crystalline SiNWs into which nickel 
atoms arc diffusing from both ends of the wire, forming 
sharp interfaces between NiSi2 "leads" and a Si "channel" 
in the center. The channel length depends on the anneal- 
ing time. Measurements of their transport characteristics 
have shown the highest on-current and on-conductance 
values recorded to date for intrinsic SiNW-FETs. It is ad- 
vantageous that the silicon channel is dopant free since 
the noise caused by structural impurities is greatly re- 
duced compared to doped nanowires. This feature is ben- 
eficial for the production of reliable SiNW-based biosen- 
sors. 

By modifying the surface of the SiNWs with 
DNAf^ii^"— enzymesi^ or other chemical compounds^ 



many kinds of biosensors such as EnFET (enzyme FETs) , 
Immuno-FETsfi^ and DNA-FETs have been demon- 
strated. Thus, SiNW-based SB-FETs are expected to 
provide a promising platform for biosensor applications. 

For a better understanding of the basic charge trans- 
port characteristics through the SiNW-based SB-FETs 
and the development of the NW-based biosensors, it is 
essential to systematically reveal the factors controlling 
the charge transport through a pristine FET device. 

Several approaches to describe the electronic trans- 
port through ID SB-FETs have been proposed. Knoch 
et al. developed quantum mechanical simulations of a 
ultrashort channel n-metail-oxide-semiconductor FETs 
(n-MOSFETs) on sihcon on insulator (SOI)^ and of 
a SB-FETs on SOI^i using real-space non-equilibrium 
Green's function formalism. Heinze et al. introduced a 
method to calculate the transport through carbon nan- 
otube (CNT) SB-FETs using the Landauer formula and 
the WKB approximation!^ Pourfath et al. extended this 
approach to cyhndrical double-gate CNT SB-FETs^ 
Michetti et al. expressed the conduction band by an ana- 
lytical formula and calculated the transmission functions 
analytically)^ Appenzeller et al. have investigated the 
transport properties of CNT-FETs and SiNW-FETs in 
detaili^ They used a modified ID Poisson equation and 
the non-equilibrium Green's function technique to calcu- 
late the transmission function with the finite difference 
method in ID. Jimenez et al. expressed the conduction 
band profile by an analytical expressio n^^'^^ and calcu- 
lated the transport properties of SB-FETs by using the 
WKB approximation, corrected with k-matching condi- 
tions for MSi2/Si(lll) and MSi2/Si(100) (with M=Ni, 
Co, and Fe) interfacesi^ 

The majority of these approaches were applied to 
highly symmetric devices such as cylindrical FETs de- 
signed for integrated circuits. In these systems, the sym- 
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metry significantly reduces the computation time. How- 
ever, to simulate realistic devices for sensor applications, 
a computationally cheap model that is applicable to ar- 
bitrary device geometries is also desired. The multiscale 
model presented in this article was developed for this 
purpose. Our method combines two approaches on differ- 
ent length scales: (1) the finite elements method (FEM) 
is used to calculate the three-dimensional (3D) electro- 
static potential across the Schottky-barrier by solving the 
Poisson equation for realistic device geometries. Then 
we extract the ID potential profile along the axis of the 
NWs and use (2) the Landauer approach combined with 
the method of non-equilibrium Green's functions^^"— to 
calculate the charge transport through the device. Our 
model is conceptually simple, computationally inexpen- 
sive, it uses only a few empirical parameters, and it is 
easily expandable. Since the method is non-atomistic the 
number of atoms in the system is irrelevant and therefore 
it is possible to simulate devices with dimensions ranging 
from a few nanometers up to some micrometers. 

In this work we model pristine SB-FETs consisting of 
SiNWs and NiSi2 NWs. We analyze the influence of the 
device geometries such as gate lengths, thickness of the 
insulators, and the gate voltages on the charge trans- 
port through the SB-FETs. Despite the simplicity of the 
model, the numerical calculations show I-V characteris- 
tics of typical conventional MOSFETs and the calculated 
saturation currents are in good agreement with the ex- 
perimental results by Weber et alt^ 

This paper is organized as follows. Section HIl presents 
theoretical framework of multiscale modeling and sec- 
tion mil shows the computational details and parame- 
ter settings. In section IIVI we present dependence of 
the electrostatic potential on device geometries, investi- 
gate transmission profiles and TV characteristics of test 
systems, and compare the numerical results with the re- 
ported experimental results. Finally, we summarize this 
paper in section |Vl emphasizing the simplicity and ver- 
satility of the multiscale modeling. 



II. THEORETICAL FRAMEWORK 

Figure[T]shows a schematic image of a SiNW-based SB- 
FET. The device consists of a SiNW working as a channel 
and NiSi2 segments working as source and drain contacts. 
This nanowire is put on the gate contact covered with 
Si02. The source, drain, and gate voltages, Vs, T^d, Vg, 
the NW diameter ^nw , as well as the length of the silicon 
nanowire channel Lc, and the thickness of the oxide layer 
tox are all indicated in figure [T] 

The band diagram of the electronic transport process 
through a SB-FET under an applied source-drain voltage 
with different gate voltages is shown in figure [5J The en- 
ergy gap $SB between the edges of the conduction band 
(CB) of the silicon channel and the Fermi energies of 
the NiSi2 contacts is called the Schottky-barrier. The 
barrier thickness is reduced with increasing gate voltage. 



allowing the electron in the source to tunnel through the 
potential energy barriers Thus, the drain current 
at the constant source-drain voltage is controlled by the 
gate voltage. The charge transport through the SB-FET 
can be divided into three stages (see figure [2]): 1- elec- 
tron (hole) injection into the conduction band (valence 
band) by tunneling through the Schottky-barrier in the 
left interface; 2. ballistic or diffusive transport of charges 
through the SiNW; 3. tunneling of charges through the 
Schottky-barrier in the right interface. 




FIG. 1: A model of a SiNW-based Schottky-barrier PET with 
realistic dimensions, as considered in this work. Vs, Vd, Vg 
are source, drain, and gate voltages, respectively. Lc is the 
length of the silicon nanowire channel, and tox is the thickness 
of the oxide layer, dnw is the diameter of the silicon nanowire. 

According to figure [2] we define the ID potential en- 
ergy profile U{x) along the axis of the NiSi2/Si/NiSi2 
nanowire as 

{A«L - *L (a: < 0) 

A^L + ^SB + C/es (0 < a; < Lc) 
ML - eV"sD - {Lc < x) 

where the position a; = is set to be at the left NiSi2/Si 
interface, ml/r is the chemical potential of the left/right 
contact, $SB is the Schottky-barrier energy, $l/r is the 
ground potential for left/right contact, Vsd is the applied 
source-drain voltage, and U-ssix) is the ID electrostatic 
potential along the axis of the Si-NW obtained from FEM 
calculations (see below) . For the Schottky-barrier energy 
$36 we use experimental values from the literature (see 
below). The ground potential $l/r is introduced in order 
to define the baseline for the transmission calculation at 
the two interfaces. The two ground potentials are cho- 
sen to span the bias window, i.e, $l = sVg^^, where 
is the maximum source-drain (bias) voltage that 
is used in a set of calculations, and $r, = 0. This is 
how the energy of the electrons in the source that are 
within bias window is set higher than the ground po- 
tential. In order to calculate the charge transmission 
at the two interfaces (see bottom of figure O we use 
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U{x) in the interval x = [— 1.5nm . . . 21.5nm] for the left 
and X = [Lc — 21.5nm . . . + l.Snm] for the right inter- 
face. In these two ranges the potential inside the SiNW 
(0 < X < Lc) is regularized by applying an energy cut- 
off J/'^"* such that U (x) = J/'^"* for energies smaller than 
t/'^"'. At the left (L) and right (R) interfaces the cutoff 
energy is U^^^ = Ml/r^^l/Rj i-S-, U{x) is never smaller 
than the ground potential $l/r (see also figure [SJ. 
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FIG. 2: An energy diagram of a SiNW-based Schottky-barrier 
FET under an applied source-drain voltage with different gate 
voltages Vg- The charge transport can be divided into 3 
processes; colierent tunneling through the Schottky-barrier 
Th/aiE) at the left/right interfaces and ballistic transport 
Tm{E) in the middle. Increasing the gate voltage leads to 
stronger bending of the conduction band (CB), which reduces 
the width of the Schottky-barrier. f&sB is the height of the 
Schottky-barrier. /il/r and ^l/r are the chemical potential 
and the ground potential for left/right contacts, respectively. 

In this study, we employed the Landauer approach in 
real-space for the calculation of the electron transmission 
through the NiSi2/Si interfaces . ^^i^^" — The Hamilton op- 
erator in ID is given by 



H = - 



2mcff dx-^ 



Uix), 



where meff , and U{x) are the effective mass of electrons 
(or holes) and the potential energy profile along the axis 
of the nanowire, respectively. Using finite differences on 
an equally spaced grid, the Hamilton matrix is given 
by23-^ 



is the potential, x„ 



where Un = U{xn) is the potential, Xn = an is the posi- 
tion, a is the grid spacing, and = It? /2'meea^ . 

The transmission functions at the left (L) and right (R) 
interfaces are obtained via the Landauer-Biittikcr formal- 



ism using the Fisher-Lee relation: 

Ti.,k{E) = Tr[G«rL/RG^r 



M 



where Tr being the trace operation and Tq {a = L, R, M) 
are the broadening functions for the contacts, given by 
Ta{E) = i\Tja{E) — Y.]^{E)]. Fm is the auxiliary broaden- 
ing function, which is required because the system is cut 
into two ID tunneling problems. The self-energies are 
defined as Yia{E) = —toexp(ika), where ka is obtained 
by inverting the band dispersion of the ID wire (linear 
chain) E{ka) = Ua + 2io(I — cos{ka)), and Ua = U{xa) 
{a = L, R, M) is the potential at the position of the 
left /right /auxiliary contactsi^"— Note that the is a 
diagonal matrix and its matrix elements are zero except 
for the position where the the left/right/auxiliary contact 
is attached. 

The retarded/advanced Green's functions G^^^ for the 
left (L) and right (R) interfaces are defined as 



G^'/^{E) = [{E±i7])I - H -J: 



L/R, 



where irj is an infinitesimal imaginary value, / is the 
identity matrix, H — Hn,m is the Hamilton matrix, and 
(a = L, R, M) are the self-energy matrices as defined 
above. 

Assuming ballistic charge transport through the silicon 
channel, i.e. Tm{E) ~ I.O, and ignoring the phase mem- 
ory, the total transmission through the device is given 
byS 

T{E) = nTn/in + Tn - TlTr). 

The assumption of ballistic transport in the core channel 
region was validated experimentally for devices having 
channel lengths less than I fimJ^ 

Finally, the current through the SB-FET is calculated 
from 

/oo 
T{E)(Fi^{E) - Fn{E)) dE. 
-oo 

The term i^L/R is the effective Fermi function /c-summcd 
over the transverse modes in the y ~ z plane, represented 

by 

Fl/r{E) = Sm,akBT/nhhn{l + exp((/XL/R - E)/kBT)), 

where S and Hl/r are the cross-sectional area of SiNWs 
and the chemical potential for left/right electrodes, re- 
spectively. Under an applied source-drain voltage VgD, 
the chemical potential for the right electrode is given by 

MR = Ml - eVsD- 

Note that the same formalism works for hole currents 
by changing the band of interest from the conduction 
band (CB) to the valence band (VB), the Fermi function 
from Fi^/^{E) to I — Fi^/^{E), effective mass of electron 
to effective mass of hole, and exchanging the sign of the 
gate voltagesi^ 



4 




FIG. 3: Electrostatic potential Ues of the SiNW-based Schottky-barrier FETs with different channel lengths being Lc — 1000 
nm in (a)/(c) and Lc = 200 nm in (b)/(d). Top: UEsix,z) in the y — plane; bottom: Ues{x) along the axes of the 
corresponding wire. In these plots the metallic NiSi2 nanowire is held at U = V. The variation of the electrostatic potential 
is strongest close to the contacts. 



III. COMPUTATIONAL DETAILS IV. RESULT AND DISCUSSION 



In this work we assumed that the effective electron 
mass of whole system for holes / electrons is equal to that 
of free electron. For all calculations we use /xl = eV, 
a=1.0 A, mcff = 1.0 X me=9. 109x10-31 kg, and T= 
300K. For the discussion of electron transport and fig- 
ure [5][H] we use tox = 300 nm, Lc = 1000 nm, cInw = 20 
nm, $SB = 0.50 eV, and V^g'jif'' = 0.5 V, unless stated dif- 
ferently. For the discussion of hole transport and figure 
[Zlilwe use iox = 300 nm, dNW = 21 nm, $sB = 0.44 eV, 
Vg"]Jf'' = 2 V. The 3D electrostatic potential UEs{x,y,z) 
for each combination of Vl, Vr, and Vg is calculated by 
solving Poisson equation with a commercial FEM soft- 
ware (COMSOL multiphysicsi^) The boundary poten- 
tials for the source, drain and gate contacts are kept con- 
stant at Vl, Vr, and Vq, respectively. The finite element 
mesh for the modeled structures is automatically gener- 
ated with controlled distribution and increased density 
in the regions close to the insulator and NiSi2/Si inter- 
faces. The ID potential profile Ubs{x) along the axis of 
the NW is extracted from UEs{x,y, z). 



A. Dependence of the electrostatic potential on 
device geometries 

As a first examination of the relationship between the 
geometry of FET devices and the efficiency of the gate 
effect, we have calculated 3D electrostatic potentials of 
FET devices. Figure EJa) and (b) depict the potential 
landscapes with a gate potential of Vg = 20 V and the 
absence of a source-drain voltage. Figure [HJc) and (d) 
show the same potential along the axes of the channels 
for long and short silicon channels, respectively. In figure 
[2a) and (c), the potential along the silicon channel drops 
strongly due to the applied gate field, whereas the gate 
field docs not penetrate efficiently in the shorter silicon 
channel in figure |3Ib) and (d). 

In addition wc have analyzed the influence of the thick- 
ness of the insulating layer between the nanowire de- 
vice and the gate on the Schottky-barrier width. The 
Schottky-barrier width Lsb of the left interface is defined 
by the position x ~ L^b > satisfying 

U{LsB)^^^L■ (1) 

In figure 21 wc show the Schottky-barrier width of the 
FET devices as a function of applied gate voltage with 
varying thickness of the gate insulators and the chan- 
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nel lengths. We can see that the gate field reduces the 
Schottky-barrier width for all devices. We find that the 
Schottky-barrier width is inversely proportional to the 
gate voltage. The reduction of the silicon channel length 
as well as an increase of the width of the oxide layer 
prevents the gate field from efficiently penetrating into 
the Schottky-contact, leaving the Schottky-barrier width 
thicker than the device with long channel or an thick ox- 
ide layer. Since the gate voltage narrows the Schottky 
barrier the enhancement of the electron tunneling by the 
gate field is larger through SB-FET devices with longer 
silicon channels and thinner gate insulators. 




I ' • ' • ^ 

4 8 12 16 20 

Gate voltage, Vq (V) 



FIG. 4: Dependence of Schottky-barrier width Lsb on the 
gate voltage Vg, the channel length Lc, and the thickness of 
the oxide layer tax- Lsb is defined in figure [1] The continu- 
ous lines are fits of the data points to Lsb = cVq^, i.e., the 
Schottky-barrier width is inversely proportional to the gate 
voltage. The gate efi'ect is most efficient for large channel 
lengths and thin oxide layers. 



B. Transmission function and I-V characteristics 

Figure [S] shows potential profiles U{x) at the left in- 
terface for different gate voltages and the corresponding 
transmission functions Ti^{E) for electron transmissions 
through these Schottky-barriers. The non-zero transmis- 
sions are due to the quantum tunneling efi'ect. Since a 
positive gate voltage causes the Schottky-barrier to get 
thinner, the electron transmission increases with increas- 
ing gate voltage. For electrons with energies above the 
Schottky-barrier height $sb, the electron transmission 
quickly approaches its maximum value 1. Therefore our 
model takes into account both tunneling currents and 
thermo-activatcd (thermionic) currents. 

After calculating the transmission functions for both 
left and right interfaces, we have calculated the I-V char- 
acteristics of SB-FET devices by integrating over the to- 
tal transmission function. Figure [6][a) shows the calcu- 
lated drain current Jgo versus source-drain voltage VgD 
for different gate voltages. The I-V curves show typical 



features of conventional FETs, i.e., a linear increase of 
current followed by current saturation for higher source- 
drain voltages. The saturation currents are strongly en- 
hanced with increasing gate fields leading to a typical 
switching behavior with a high on/off-current ratio. We 
evaluated this behavior by calculating the drain current 
at a fixed source-drain voltage for several gate voltages. 
Figure [Bib) presents the saturated drain current versus 
the gate voltage for a fixed source-drain voltage Vsd = 0.5 
V. The drain current increases exponentially with the in- 
crease of the gate field in the lower gate voltages. This is 
due to the rapid reduction of the Schottky-barrier width 
with the increase of gate fields as shown in the figure 
m The saturation of the drain current in higher gate 
field can also be explained from the figure |4] since the 
Schottky-barrier width does not decrease so much in the 
higher gate field. The huge on/off-currcnt ratio implies 
that the device is promising for logic operations. The 
flatness of the current in the saturation region in figure 
inia) is also beneficial for the logic operation since the 
current retains constant with respect to fluctuation of 
the source-drain voltage. 



C. Comparison with experiment 

Many measurements of SB-FETs consisting of 
SiNWs have exhibited unipolar p-type transfer 
characteristics,— "—>2^ thus the current is transported 
by holes in these systems. As discussed previously, our 
model can also be applied to hole transport systems. In 
order to check the validity of the model, we have applied 
our model to SiNW-based SB-FETs and compared the 
I-V characteristics and on/off-current ratios with the 
experimental results measured by Weber et al^ 

The numerical parameters are set as follows. The 
Schottky-barrier height between Ni-silicide and bulk Si 
ranges from 0.35 eV to 0.47 eV for holes.— Thus, we set 
the Schottky-barrier height of the NiSi2/Si interface for 
holes as <f>sB= 0.44 cV. The diameter of the SiNWs is 
set to c?NW = 21 nm, identical to the mean diameter of 
the literaturCfi^ The oxide thickness on the gate is set to 
tox =300 nm. 

Figure[7]presents the calculated drain currents /gD ver- 
sus source-drain voltage T^sd for (a) long (Lc = 1000 nm) 
and (b) short gate length (Lc = 200 nm) with experimen- 
tal data provided by Weber et alr^ The calculated I-V 
curves for both systems exhibit similar features as the ex- 
perimental results having both a linear increase and sat- 
uration at higher source-drain voltages. The drain cur- 
rent through the shorter channel is smaller than the one 
through the longer channel since the gate field does not 
optimally penetrate the NiSi2/Si interfaces. Although 
the slopes of the calculated TV curves for low source- 
drain voltages significantly exceed experimental ones, the 
amount of the saturated drain current shows a good 
agreement with the experiment except for the currents 
with high gate field. 
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FIG. 5; (a) The potential profile 17(2;) across the central axis of the left NiSi2/Si interface for different gate voltages Vg, and 
(b) the corresponding transmission functions Tl{E). Here the Schottky-barrier height $sb = 0.50 eV. 
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FIG. 6: Electron transport characteristics of the SB-FET device, (a) Drain current Isd versus source-drain voltage Vsd for 
different gate voltages Vg- The drain current saturates at higher bias voltages. The saturated current is significantly enhanced 
with higher gate voltages, (b) Drain current 7sd versus gate voltage Vg for a fixed source-drain voltage (Vsd ~ 0.5 V). We 
obtain the typical behavior of a conventional FET device. 



Figure [5] presents the saturated drain current versus 
the gate voltage of SB-FET devices with long and short 
channels. The saturated drain currents increase signifi- 
cantly with the increase of the gate field since the gate 
field reduces Schottky-barrier width as shown figure SI 
The on-currents in the short and long systems are of the 
same order of magnitude. On the other hand, the off- 
currents for small gate voltages are significantly smaller 
for the short channel. Although the drain current in the 
lower gate field is underestimated, the numerical results 
show a qualitative agreement with experiment. Generally 
we observe that the drain currents of the SB-FET with 
the long channel are higher than for the short channel 
because the gate effect is more efficient (compare figure 
He), (d), figure a and figure [Ua), (b)). 

The discrepancies between our model and the experi- 
ment come from the simplicity of our model. First, we 
did not take into account the sub-bands^li^^i^ but ex- 
pressed the transport by using a single band model. The 
inclusion of the sub-bands will lower the drain current 
for small source-drain voltages since the onset of cur- 
rent flow through the channel will be shifted to higher 



source-drain voltages due to the difference of energies be- 
tween sub-bands. In addition, the reconstruction of the 
exact potential profile at the Schottky-contac t^°i^^ due to 
space-charge effects and the dependence of the potential 
on the lateral position in the wire are not considered in 
this study. Furthermore, the effect of capacitance at the 
left and right interfaces and at the gate contact are not 
includedJ^, The addition of a native Si02 layer on the 
surface of the nanowire and the corresponding change of 
the electronic structures between the core and surface of 
the nanowire will also have an impact on the transport 
properties!^ In order to capture the underlying physics 
of SB-FET devices, an advanced model incorporating 
these effects will be needed. Further studies including 
these effects will give helpful information for the future 
development of nanowire-based FET devices. 



V. CONCLUSION 

In summary, nanowire-based FETs have been focus- 
ing attention as a promising platform for sensor applica- 
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FIG. 7: I-V characteristics of two FET devices for (a) a long 
and (b) a short channel length at different gate voltages. The 
calculated currents are shown with dots. The experimental 
results given by Weber et al;^ are shown as continuous lines 
for comparison. Although the theoretical currents reach their 
saturations at lower bias, the saturated currents show a good 
agreement with experiments. Note that Ids = — ^sd and 
Vbs = — VsD. 
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FIG. 8: Drain current Ids versus gate voltage Vg character- 
istics of SiNW-based SB-FET devices with long (in red) and 
short (in green) silicon channels. The source-drain voltage 
is fixed to be Vos ~ 0.5 V. For comparison, the experimen- 
tal resuhs of Weber et alM for a long (1300 nm) SiNW are 
shown (in blue). Numerical results show a qualitative agree- 
ment with experiment. The drain currents of the SB-FET 
with the long channel are higher than for the short channel 
because the gate effect is more efficient. 



Our approach is suitable for one-dimensional Schottky- 
barrier field-effect transistors of arbitrary device geome- 
try. The model can be improved by taking into account 
for example insulating layers surrounding the nanowires, 
the deformation of the electrostatic potential at the 
Schottky-contacts, or accurate electronic band disper- 
sions. The results obtained from our model will serve as 
guidcUncs for the understanding of ongoing experimental 
work and will help in the design of device architectures 
for future applications. 



tions due to their high sensitivity. In order to understand 
the origin of the physical behavior in the nanowire-based 
FET devices, we have developed a multiscale model com- 
bining classical FEM and the Landauer approach. Our 
model is conceptually simple, computationally inexpen- 
sive, it uses only a few empirical parameters, and it is pos- 
sible to simulate devices with dimensions ranging from a 
few nanometers up to some micrometers. We have ap- 
plied this model to Schottky-barrier FETs consisting of 
SiNWs working as channel and NiSi2 NW working as 
source and drain contacts. Our calculated I-V character- 
istics showed the typical behavior of conventional FET 
devices, having a linear increase of current followed by 
a saturation. The saturated drain currents increase with 
increasing gate field due to the reduction of the Schottky- 
barrier width. Despite the simplicity of the model, the 
calculations showed a good agreement with experiments 
and the model correctly reproduces the dependence of the 
saturated drain current on the gate field and the device 
geometry. 
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